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ACHILL SCHURMANN 

Abstract. A positive definite quadratic form is called perfect, if it is 
uniquely determined by its arithmetical minimum and the integral vec- 
tors attaining it. In this self-contained survey we explain how to enu- 
merate perfect forms in d variables up to arithmetical equivalence and 
scaling. We put an emphasis on practical issues concerning computer 
assisted enumerations. For the necessary theory of Voronoi we provide 
complete proofs based on Ryshkov polyhedra. This allows a very natu- 
ral generalization to T-perfect forms, which are perfect with respect to a 
linear subspace T in the space of quadratic forms. Important examples 
include Gaussian, Eisenstein and Hurwitz quaternionic perfect forms, 
for which we present new classification results in dimensions 8, 10 and 
12. 



1. Introduction 

In this paper we are concerned with perfect forms, which are real positive 
definite quadratic forms 

d 

(1) Q[x]= qijXiXj 

in d variables x = (x\, . . . E determined uniquely by their arith- 
metical minimum 

(2) A(Q)= nun ON 

x<=Z d \{0} 

and its representations 

(3) Min Q = {x eZ d : Q[x] = A(Q)}. 

The study of perfect forms goes back to the work of Korkin and Zolotarev 
|KZ77j . They observed that perfection is necessary for positive definite 
quadratic forms in order to give a local maximum of the Hermite invariant 

Such forms are called extreme. 
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As briefly reviewed in Section [21 finding the global maximum of the Her- 
mite invariant, or equivalently the densest lattice sphere packing is a widely 
studied problem. In this article we describe the only known algorithmic 
solution of this problem which works in principle in every dimension. It is 
based on the classification respectively enumeration of perfect forms. We 
refer to |RB79j . [CS99] . [Mar03j . [GruC>7] and [Sch08| for further reading. 

Based on perfect forms, Voronoi [Vor07] developed a polyhedral reduction 
theory, which was later found to have several applications in other contexts. 
It has for example been used for compactification of moduli spaces (cf. for 
example |AMRT75| . |McC98j . |She06 j). for computing the cohomology of 
GI_d(Z) and of congruence subgroups, as well as for computing algebraic 
if -groups K<i(Z) for small d and up to small torsion (cf. |Sou99j . [EGS02J 
and the appendix in [Ste07| ). A basic task in these computations is the 
enumeration of perfect forms. In some of the applications it is also necessary 
to understand more of the structure of the Ryshkov polyhedron (to be defined 
in Section [3]) whose vertices are perfect forms. 

In this article we explain Voronoi's theory based on the Ryshkov polyhe- 
dron. We provide complete proofs for all of its required properties. We think 
that this view is more accessible than the usual dual viewpoint, originally 
taken by Voronoi and by most other authors subsequently. Voronoi's algo- 
rithm can be simply described as a traversal search on the graph consisting 
of vertices and edges of the Ryshkov polyhedron. This viewpoint allows 
in particular a very simple and direct generalization to so called T-perfect 
forms: Intersecting a linear subspace T with the Ryshkov polyhedron yields 
a lower dimensional Ryshkov polyhedron whose vertices are T-perfect forms. 
Voronoi's theory immediately generalizes. 

The article is organized as follows. In Section [2] we review some necessary 
background and notations. In Section [3] we define the Ryshkov polyhedron 
and prove that it is "locally finite". This yields the grounds for Voronoi's 
algorithm to be described in Section HI Here we put special emphasis on 
practical issues related to running Voronoi's algorithm on a computer. In 
Section [5] we briefly explain how to determine extreme forms. Section [6] 
contains some informations on automorphism groups and their computation 
and in Section [7] we explain the "T-theory" , when restricting to a linear 
subspace T. As examples of linear subspaces that contain forms invariant 
with respect to a finite group of automorphisms, we consider in Section [8] 
forms with a Gaussian, Eisenstein or Hurwitz quaternionic structure. We 
obtain several new classification results. 



2. Background on positive definite quadratic forms 

In this section we review - basically from scratch - some of the historical 
background and notations used in the remaining of the article. The reader 
familiar with most of this background may simply skip this section. 
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We consider real quadratic forms in d variables as in ([I]), hence with 
coefficients qij G R. By assuming qij = qji without loss of generality, we 
simply identify the quadratic form Q with the real symmetric matrix Q = 
(Qij)i,j=i,...,d- The space of all real quadratic forms in d variables is identified 
with the space 

S d = {Q G R dxd : Q* = Q} 

of real symmetric d x d matrices. Using matrix notation we have Q[x] = 
x l Qx. Endowed with the inner product 

d 

(Q, Of) = Yl = tiac <Q ■ <9'), 

S d becomes a C^ 1 ) -dimensional Euclidean space. 

Two quadratic forms Q,Q f G 5 rf are called arithmetically (or integrally) 
equivalent, if there exists a matrix f7 in the group 

Gl_d(Z) = {[/ G Z dxd : | det 17| = 1} 

such that 

Q' = U l QU. 

Note that = Q'\L d \ for arithmetical equivalent Q and Q', but the 

opposite may not hold. 

A quadratic form Q G S d is positive definite, if > for all x G 

R \ {0}. The set of all positive definite quadratic forms (PQFs from now 
on) is denoted by 5> . It is not hard to see that S d is an open (full 
dimensional) convex cone in S d with apex 0. In particular for Q G S d , 
the open ray {XQ : A > 0} is contained in S d as well. Only for PQFs the 
arithmetical minimum defined in ((2J) is greater than 0. 

A PQF Q defines a real valued strictly convex function on R d and for 
A > 

(5) E(Q,\) = {x G R d : Q[x] < A} 

is a non-empty ellipsoid with center 0, providing a geometric interpretation 
of a PQF. The arithmetical minimum is the smallest number A > for which 
the ellipsoid E(Q, A) contains an integral point aside of 0. The integral 
points x in MinQ (see (J3j) ) lie on the boundary of the ellipsoid E(Q, A(Q)). 

Hermite, who initiated the systematic arithmetic study of quadratic forms 
in d variables found in particular an upper bound of the arithmetical mini- 
mum in terms of the determinant det Q of Q: 

Theorem 1 (Hermite, |Her50| ). 

/ 4 n (d-l)/2 

A(Q) < (detQ) 1 ^- - forallQeS d . 
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d 


lattice 


5d 


'rid 


author (s) 


n 

L 


A„ 
"2 


u.yuoy . . . 






3 


A 3 = D 3 


0.7404 . . . 


2 l/3 


Gauss, 1840, [Gau40j 


4 


D 4 


0.6168... 


4l/4 


Korkin & Zolotarev, 1877, |KZ77j 


5 


D 5 


0.4652 . . . 


8 V5 


Korkin & Zolotarev, 1877, |KZ77j 


6 


E 6 


0.3729... 


(f) V6 


Blichfeldt, 1935, |Bli35j 


7 


E 7 


0.2953... 


64 l/7 


Blichfeldt, 1935, [Bh35] 


8 


E 8 


0.2536... 


2 


Blichfeldt, 1935, |Bli35j 


24 


A 24 


0.0019... 


4 


Cohn & Kumar, 2004, |CK09j 



Table 1. Known values of Hermite's constant. 



Hermite's theorem implies in particular the existence of Hermite's con- 
stant 

Hermite's constant and generalizations have been extensively studied, e.g. 
in the context of algebraic number theory and differential geometry. We 
refer to |Bav97j . |Sch98j . |Cou01j and |Wat04j for further reading. 

The following lattice sphere packing interpretation is due to Minkowski: 
Using a Cholesky decomposition Q = A f A of a PQF Q, with A € GL^(IR), 
the set L = AZ d is a (point) lattice, that is, a discrete subgroup of M. d . The 
column vectors of the matrix A are referred to as a basis of L. The maximum 
radius of non overlapping solid spheres around lattice points of L is 

the so called packing radius of L. Denoting the solid unit sphere by B d , the 
sphere packing density 6(L) of a lattice L is defined as the portion of space 
covered by solid spheres of radius A(L), hence 

N vol(A(L)£ d ) X(L) d volB d 
d(L) = 



det L det L 

Note that J is invariant with respect to isometries and scalings of the lattice 
L. The supremum of possible lattice packing densities 5d is, up to a constant 
factor, equal to a power of Hermite's constant. Table [T] lists the dimensions 
in which 5d respectively Hermite's constant fid is known. 

The lattices for d > 2, for d > 3 and for d = 6, 7, 8 are the 
so-called root lattices. One of the most fascinating objects is the Leech 
Lattice A24 in 24 dimensions. Definitions and plenty of further information 
on these fascinating lattices can be found in jCS99j . jMar03l and the online 
database INSI. 
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Minkowski noticed [Min91| that the trivial bound 

(7) 5(L) < 1, 

which is an immediate consequence of the sphere packing interpretation, 
tremendously improves the upper bound for the arithmetical minimum in 
Hermite's Theorem. In fact, ([7|) is equivalent to 

(8) MQ) < (det Q)V,._±_. 

showing that the exponential constant on the right in Theorem Q] can be 
replaced by a constant which grows roughly linear with d. 

This trivial, but significant improvement lead Minkowski to a powerful 
fundamental principle. The ellipsoid E(Q, rg), with tq being the right hand 
side in ([5]), has volume 

vol E(Q,r Q ) = vol{^A- l B d ) = r^ 2 (det Q)~ 1/2 vol B d = 2 d . 

Minkowski discovered that not only ellipsoids of volume 2 d contain a non- 
zero integral point, but also all other centrally symmetric convex bodies 
(non-empty, compact convex sets). 

Theorem 2 (Minkowski's Convex Body Theorem). Any centrally symmet- 
ric convex body in M. d of volume 2 d contains a non-zero integral point. 



3. Ryshkov POLYHEDRA 

Since the Hermite invariant is invariant with respect to scaling, a natural 
approach of maximizing it is to consider all forms with a fixed arithmetical 
minimum, say 1, and minimize the determinant among them. We may even 
relax the condition on the arithmetical minimum and only require that it is 
at least 1. In other words, we have 

H d = l/inf(detQ) 1/d , 

where 

(9) K={QeS d :\(Q)>l}. 

We refer to 1Z as Ryshkov polyhedron, as it was Ryshkov |Rys70| who no- 
ticed that this view on Hermite's constant allows a simplified description of 
Voronoi's theory. 

Because of the fundamental identity 

Q[x) = (Q,xx f ), 

quadratic forms Q G S d attaining a fixed value on a given x G W d \ {0} lie all 
in a hyperplane {affine subspace of codimension 1). Thus Ryshkov polyhedra 
1Z are intersections of infinitely many half spaces: 

(10) K = {Q G S d : (Q, xx l ) > A for all x G Z d \ {0}}. 
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We show below that are 1Z is "locally like a polyhedron" . Its vertices are 
precisely the perfect forms with arithmetical minimum 1. 

Background on polyhedra. Before we give the precise statement, and for 
later purposes, we need some basic notions from the theory of polyhedra. 
As general references for further reading we recommend the books [MS71| . 
[Shr86], [Zic97j, |Grii02]. A convex polyhedron V C E in a Euclidean space 
E with inner product (•,•) (e.g. E = S d ) can be defined by a finite set of 
linear inequalities (TC- description) 

V = {x £ E : (aj, x) > bi, % = 1, . . . , to}, 

with en £ E and i>, £ R for i = 1, . . . ,m. If the number of inequalities to 
in the description is minimum, we say it is non-redundant. The dimension 
dim"P of V is the dimension of the smallest affine subspace containing it. 
Under the assumption that V is full-dimensional every inequality i of a non- 
redundant description defines a facet {x £ V : (ai,x) = bi} of V, which is 
a (d — l)-dimensional convex polyhedron contained in the boundary of V . 
More generally, an intersection of a hyperplane with the boundary of V is 
called a face of V, if V is contained in one of the two halfspaces bounded by 
the hyperplane. The faces are polyhedra themselves; faces of dimension 
and dimension 1 are called vertices and edges. 

By the Farkas-Minkowski-Weyl Theorem (see e.g. [Shr86l Corollary 7.1a]), 
V can also be described by a finite set of generators (V '-description): 

V = couv{ui,...,Ufe} + cone{ufe + i,...,u n } 

n k 

= {X^:A;>0,^Ai = l} 

i=l i=i 

where v j £ E for i = 1, . . . ,n. Here convM denotes the convex hull and 
coneM the conic hull of a set M. If the number of generators is minimum, 
the description is again called non-redundant. In the non-redundant case, 
the generators V{, i = 1, . . . , k, are called vertices and R>o"Vj, i = k+1, . . . , n, 
are the extreme rays of "P. In case V is bounded we have n = k and we speak 
of a convex poly tope. 

There exist several different approaches and corresponding software for 
the fundamental task of converting ^-descriptions of polyhedra into V- 
descriptions and vice versa (see for example cdd [SoCdj and Irs [SoLrj ). 

Locally finite polyhedra. We say that an intersection of infinitely many 
halfspaces, V = Cl^iH^, is a locally finite polyhedron, if the intersection 
with an arbitrary polytope is a polytope. So, locally V "looks like a poly- 
tope" . 

Theorem 3. For d > 1, the Ryshkov polyhedron 1Z (see (flOl) ) is a locally 
finite polyhedron. 
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Proof. By applying Minkowski's convex body Theorem [21 we show below 
that 

(11) TZn{Q G S d : trace Q < C] 

is a polytope (possibly the empty set) for every constant C. This proves the 
theorem, since 

trace Q = (Q, id^) < C 

determines a halfspace containing a bounded section of 5> - 

The sets (jlip are polytopes if the set of all x G Z d \ {0} with Q[x] = 1 
(or Q[x] < 1) for some forms Q in (jlip is finite. We show below that the 
absolute value of coordinates 

m = max \x{\ 

i=l,...,d 

of x with this property is bounded. 

Let Q be a PQF in p}. Then the ellipsoid £(Q, 1) = {x G M d : < 1} 
does not contain any point of Z d \ {0} in its interior. So in particular 
vol E(Q, 1) < 2 d by Minkowski's convex body theorem. Since 

1 < Q[ei] < (trace Q) - ^ Q[ej] < C - (d - 1), 

we know that i?((5, 1) contains the cross polytope 

(12) C"-conv{± ei : i = l,...,rf} 
with 

C" = (C- (d- 1)) _1/2 . 

For x with < 1 consider the polytope defined as the convex hull of ±x 
and the cross polytope (|12p . It is contained in E(Q, 1). On the other hand, 
this polytope contains the convex hull P of ±x and the (d — l)-dimensional 
cross polytope 

C ■ conv{±ej : i = 1, . . . , d, i ^ j}, 

where j G {1, . . . ,n} is chosen such that \xj\ attains m. Thus setting C" 
to be the (d — l)-dimensional volume of latter (d — l)-dimensional cross 
polytope we get 

m ■ -C" = volP < vol (conv{±x, Rgjft) < vol E(Q, 1) < 2 d . 
d 

Hence we obtain the desired bound on m (depending only on d). □ 

One consequence of the Theorem is the fact that Hermite's constant can 
only be attained by perfect forms, which was first observed by Korkin and 
Zolotarev in [KZ77]. This follows immediately from the following Theorem. 

Theorem 4 (Minkowski [Mm05]). (detQ) 1 ^ is a strictly concave function 
on S d Q . 
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For a proof see for example [GL87j . Note, that in contrast to (det Q) 1 ^, 
the function det Q is not a concave function on S d (cf. |Nel74] ) . However 
Minkowski's theorem implies that the set 

(13) {Q G S d : det Q > D} 

is strictly convex for D > 0. 

Finiteness up to equivalence. The operation of GL^(Z) on S d leaves 
X(Q), MinQ and also 1Z invariant. GL^(Z) acts on the sets of faces of a 
given dimension, thus in particular on the sets of vertices, edges and facets 
of 1Z. The following theorem shows that the Ryshkov polyhedron 1Z contains 
only finitely many arithmetically inequivalent vertices. By Theorem [J] this 
implies in particular that Ti^ is actually attained, namely by some perfect 
forms. 

Theorem 5 (Voronoi 1907). Up to arithmetical equivalence and scaling 
there exist only finitely many perfect forms in a given dimension d > 1. 

Proof. In the proof of Theorem [3] we showed that the set (|lip of PQFs 
Q with X(Q) > 1 and trace Q < C is a polytope, hence has only finitely 
many vertices. Therefore it suffices to show that every perfect PQF Q 
with X(Q) = 1 (a vertex of the Ryshkov polyhedron 1Z) is arithmetically 
equivalent to a form with trace smaller than some constant depending only 
on the dimension d. 

By Hermite's Theorem [T] we find an equivalent PQF Q' with 

d /4 xd(d-l)/2 

(14) u - det ^- 

1=1 ^ ' 

The determinant det Q' = det Q can be bounded by 1 because of Hadamard's 
inequality showing 

(15) detg< Q[ax]---Q[a d ] 

for Q £ 5> and linearly independent a\,...,ad € Z d . Latter applies in 
particular to linearly independent vectors in MinQ, respectively MinQ'. 
The existence of d linear independent vectors in Min Q for a perfect form Q 
follows from the observation that the rank-1 forms xx l with x S MinQ have 
to span S d , since they uniquely determine Q through the linear equations 
(Q,xx*) = X(Q). If however MinQ does not span M. d then these rank-1 
forms can maximally span a ( 2 ) -dimensional subspace of S d . 
Because of q' u > 1 we find 

Qkk < II q 'a - ( i 

i=i \' 

From this we obtain the desired upper bound for the trace of Q'\ 

d /ylN d(d-i)/2 



/4\ ' 

trace Q' = Y^l'kk < d (j^J 
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□ 

4. VORONOl'S ALGORITHM 

The vertices (perfect PQFs) and edges of TZ form the (abstract) Voronoi 
graph in dimension d. Two vertices, respectively perfect PQFs Q and Q' 
are connected by an edge if the line segment conv{Q, Q'} is an edge of 1Z. 
In this case we say that Q and Q' are contiguous perfect forms (or Voronoi 
neighbors). By Theorem [5l for given d, there are only finitely many vertices 
(and edges) of the Voronoi graph up to arithmetical equivalence. Therefore, 
one can enumerate perfect PQFs (up to arithmetical equivalence and scaling) 
by a graph traversal algorithm, which is known as Voronoi 's algorithm (see 
Algorithm [1]) . 



Input: Dimension d. 

Output: A complete list of inequivalent perfect forms in 5> . 
Start with a perfect form Q. 

1 . Compute Min Q and describing inequalities of polyhedral cone 
(16) V{Q) = {Q 1 G S d : Q'[x] > for all x G MinQ} 

2. Enumerate extreme rays R\, . . . ,Rj. of the cone ^(Q) 

3. Determine contiguous perfect forms Qi = Q + aRi, i = 1, . . . , k 

4. Test if Qi is arithmetically equivalent to a known form 

5. Repeat steps 1.-4. for new perfect forms 

Algorithm 1. Voronoi's algorithm. 

As an initial perfect form we may for example choose Voronoi's first 
perfect form, which is associated to the root lattice A^. For example take 
Q Ad = (Qi,j )i<i,j<d with q iti = 2, = <?j_ M = -1 and q itj = otherwise 

(see Section 6.1] or |Mar03l Section 4.2]). 

One key ingredient, not only for step 1., is the computation of represen- 
tations of the arithmetical minimum. For it we may use the Algorithm of 
Fincke and Pohst (cf. |Coh93] ): Given a PQF Q, it allows to compute all 
x G 7h d with Q[x] < C for some constant C > 0. For C = min, = i viii rf qa 
a non-zero integral vector x with Q[x] < C exists, hence in particular 
X(Q) < C. The Fincke and Pohst algorithm makes use of the Lagrange 
expansion of Q, given by 

d I d 
(17) Q[x] =y^ Aj \xj- y~] aijXj 

i=i y i=i+i 

with unique positive outer coefficients A4 and inner coefficients aij G M, for 
i = l,...,d and j = i + 1, . . . , d. By it, it is possible to restrict the search 
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to integral vectors x with 



d 





OtijXj 



< 



1=1 



for i = d, . . . , 1. Here, the bound on the coordinate Xi depends on fixed val- 
ues of Xj+i, . . . , Xd, for which we have only finitely many possible choices. Im- 
plementations are provided in computer algebra systems like Magma [SoMaJ 
or GAP [SoGa| (see also shvec by Vallentin SoSh]). 

For step 2., observe that the homogeneous cone (|16p is a translate of the 
support cone 



of Q at TZ. Having its W-description (by linear inequalities) we can transform 
it to its V-description and obtain its extreme rays. The extreme rays R 
provided by Q through (fT6j) are easily seen to be indefinite quadratic forms 
(see [Mar03]). 

In step 3., the contiguous perfect forms (Voronoi neighbors) of Q are 
of the form Q + pR, where p is the smallest positive number such that 
A(Q + pR) = A and Min(Q + pR) ^ MinQ. It is possible to determine p, 
for example with Algorithm [2j 



Input: A perfect form Q G 5> and an extreme ray R of (|16p 
Output: p > with X(Q + pR) = \{Q) and Min(Q + pR) % MinQ. 

(l,u) <- (0,1) 

while Q + uR(£ S* or X(Q + uR) = \{Q) do 

if Q + uR then u <- (I + u)/2 

else (/, u) <— (u, 2u) 

end if 
end while 

while Min(Q + IR) C MinQ do 



{Q' G 5 d : Q'[x] > Q[s] for all x G MinQ} 




if A(Q + 712) > A(Q) then I <- 7 



else 



u «- min{(A(Q) - Q[v])/R[v] : v G Min(Q + 7 12), R[v] < 0} U {7} 



end if 



end while 



P 



Algorithm 2. Determination of Voronoi neighbors. 



In phase I (first while loop), the procedure determines lower and upper 
bounds I and u for the desired value p, such that Q + IR, Q + uR G 5> 
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r] 

fjj 


-f-f- Tiprf fnrmQ 

-f-f- pen. J-VJI lllO 


-f-f- fi"v"i" fmrnnc! 

-H— C-A. \j • L\JL lllo 


£1 1 1 1 h OT* 1 t; 1 


2 


1 


1 


Lagrange, 1773, Lag 73 


3 


1 


1 


Gaufi, 1840, [Gau40] 


4 


2 


2 


Korkin & Zolotarev, 1877, |KZ77j 


5 


3 


3 


Korkin & Zolotarev, 1877, |KZ77j 


6 


7 


6 


Barnes, 1957, [Bar57j 


7 


33 


30 


Jaquet-Chiffelle, 1993, |Jaq93| 


8 


10916 


2408 


Dutour Sikiric, Schiirmann & Vallentin 


9 


> 500000 




2005. |SV05|.[DSV07|. cf. |Rie06j 



Table 2. Known numbers of perfect and extreme forms. 



with X(Q + IR) = A and X(Q + uR) < A. In phase II, the value of p is deter- 
mined. Note that replacing the assignment of u by the simpler assignment 
u <— 7 corresponds to a binary search coming at least arbitrarily close to p. 
However, it may never reach the exact value. 

For step 4. observe, that based on an algorithm to compute short vectors 
(for example the one by Fincke-Pohst described above), it is possible to test 
algorithmically if two PQFs Q and Q' are arithmetically equivalent. That 
is, because the existence of U £ GL^(Z) with Q' = U l QU implies 

Qu = Q'[<H] = Q[ui\- 

Hence for the i-th column Ui of U we have only finitely many choices. This 
idea, but more sophisticated, is implemented in isom by Plesken and Sou- 
vignier [PS97], which is also part of Magma [SoMa] and Carat [SoCaj. Note 
that isometry tests for perfect forms can be simplified, because it suffices to 
find a U G Gl d (Z) with C/MinQ' = MinQ. 

Using the described software tools it is possible to verify the results of 
Table [2] below on any standard PC up to dimension 6. Note however, that 
this computation was already done without a computer by Barnes |Bar57| . 
In dimension 7 and beyond the explained procedure has a seemingly insu- 
perable "bottleneck": The enumeration of extreme rays for support cones 
with many facets, respectively for perfect forms with large sets MinQ- 

There have been several attempts of using computers to (try to) enumer- 
ate perfect forms. Larmouth [Lar71] was the first who implemented it and 
was able to verify the result of Barnes [Bar57] up to dimension 6. Also, 
Stacey [Sta75j and Conway and Sloane [CS88] used computer assistance for 
their attempts to classify the perfect forms in dimension 7. Exploiting sym- 
metries, Jaquet-Chiffelle |Jaq93| was able to enumerate all perfect forms in 
dimension 7. Recently, together with Mathieu Dutour Sikiric and Frank 
Vallentin we were able to finish the classification in dimension 8 (see [SV05| 
and [DSV07] ). 
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5. EUTAXY AND EXREMALITY 

Not every perfect form is extreme, hence gives a local maximum of the 
Hermite invariant, as shown in Table [2] from dimension 6 onwards. 

In order to characterize extreme forms the notion of eutaxy is used: A 
PQF Q is called eutactic, if its inverse Q^ 1 is contained in the (relative) 
interior relint V(Q) of its Voronoi domain 

V(Q) = cone-jW : x G MinQ}. 

Note that the Voronoi domain is full-dimensional if and only if Q is perfect. 
Note also that the rank-1 forms xx t give inequalities (Q, xx t ) > 1 defining 
the Ryshkov polyhedron and by this the Voronoi domain of Q is equal to 
the normal cone 

(18) {N £S d : (N, Q/X(Q)) < (N, Q') for all Q' G 11} 
of U at Q/X(Q). 

Algebraically the eutaxy condition Q~ l G relint V(Q) is equivalent to the 
existence of positive a x with 

(19) Q~ l = ot x xx l . 

xGMinQ 

Computationally, eutaxy of Q can be tested by solving the linear program 

(20) maxa m i n s.t. a x > a m i n and (Tl9|) holds. 

The form Q is eutactic, if and only if the maximum is greater 0. 

Voronoi [Vor07] showed that perfectness, together with eutaxy implies 
extremality and vice versa. (Eutaxy alone does not suffice for extremality.) 
By solving the linear program (|20p for perfect forms a list of extreme forms 
can be obtained. This was done by Riener [Rie06] for the 8 dimensional 
perfect forms, showing that only 2408 of them are extreme (see Table [2]). 

Geometrically the characterization of extreme forms by Voronoi can easily 
be seen from the identity 

(21) grad det Q = (det Q)Q~ l 

for the gradient of det Q. By it, the tangent hyperplane T in Q of the smooth 
determinant-det Q -surface 

S = {Q' G S* : det Q' = det Q} 

is given by 

T = {Q 1 G S d : (Q-\ Qf) = (Q-\Q)}. 

Or in other words, Q _1 is a normal vector of the tangent plane T of S at 
Q. By Theorem |4] the surface S is contained in the halfspace 

(22) {Q'£S d : {Q-\Q'-Q) >0}, 

with Q being the unique intersection point of S and T. 

As a consequence, a perfect form Q with X(Q) = 1 attains a local mini- 
mum of det Q (hence is extreme) if and only if the halfspace (122j) contains 
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the Ryshkov polyhedron 1Z, and its boundary meets 1Z only in Q. This is 
easily seen to be equivalent to the condition that the normal cone (Voronoi 
domain) V(Q) of 1Z at Q contains Q^ 1 in its interior. 

6. Automorphism groups 

The recent enumeration success in dimension 8 was previously not pos- 
sible, because the computation of extreme rays was in particular difficult 
for the support cones associated to the highly symmetric forms associated 
to the root lattices E7 and Eg. Note that the enumeration of extreme rays 
is a known difficulty in many problems, for example in combinatorial op- 
timization. Martinet stated that "it seems plainly impossible to classify 
8-dimensional perfect lattices" (see [Mar 031 p.218]). However, it is possible 
to overcome these difficulties to some extend by exploiting symmetries in 
the computation. For a survey on such symmetries exploiting techniques we 
refer to |BDS09j . 

In general the automorphism group (or symmetry group) of a quadratic 
form Q G S d , is defined by 

AutQ = {U G GL d (Z) : U l QU = Q}. 

As in the case of arithmetical equivalence, we can determine Aut Q, based on 
the knowledge of all vectors u G Z d with Q [u] = qu for some i G {1, . . . , d}. 
Again, Magma [SoMaj, based on an implementation of Plesken and Souvignier 
(also available in Carat [SoCaJ), provides a function for this task. 

For Q G 5> with X(Q) = 1, the support cone V(Q) at Q of the Ryshkov 
polyhedron 1Z (see ([16]) ) and its dual, the Voronoi domain V(Q), inherit 
every symmetry of Q. That is, for all U G Aut Q we have 

XJ l V{Q)U = V(Q) and U l V(Q)U = V(Q). 

The automorphism group of a PQF Q is always finite. On the other 
hand, for every finite subgroup G of GL^(Z), there exists a PQF Q with 
G C Aut Q. For example, given an arbitrary Q' G S^ , the PQF 

Q=Y1 ut Q' u 

UeG 

is invariant with respect to G, hence satisfies G C Aut Q. 

For a finite group G C GL^(Z), the space of invariant quadratic forms 

(23) T G = {Q G S d : U l QU = Q for all U G g| 

is a linear subspace of S d ; Tq n 5> is called Bravais manifold of G. 

7. T-PERFECT FORMS 

Since the enumeration of all perfect forms becomes practically impossible 
in higher dimensions (due to the complexity of the Ryshkov polyhedron 1Z), 
it is natural to restrict classifications to certain Bravais manifolds. This is in 
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particular motivated by the fact that all forms known to attain the Hermite 
constant have large symmetry groups. 

Within Tq we are lead to the theory of G-perfect forms of Berge, Martinet 
and Sigrist jBMS92j , It generalizes to a theory of T -perfect forms, where 
T C S d is some linear subspace (see [Mar03| ). Suitable linear subspaces 
T allow systematic treatments of important classes of forms. Examples 
are Eisenstein, Gaussian and Hurwitz quaternionic forms as explained in 
Section [8l For further informations on classes as cyclotomic forms or forms 
having a fixed section we refer to |SigOO| and [Mar03j . 



Our viewpoint developed in this article (based on Ryshkov polyhedra) 
allows a straightforward description of the "T-theory" . Given a linear sub- 
space T C S d we simply consider the intersection 

(24) TZf)T. 

It is again a locally finite polyhedron which we call a Ryshkov polyhedron 
too. Its vertices are called T -perfect forms. In case T = Tq, where G is a 
finite subgroup we speak of G-perfect forms. One should be aware that in 
general, T-perfectness does not imply perfectness. 

We have to modify the notion of equivalence. Two PQFs Q and Q' are 
called T-equivalent if there exists a U S GLrf(Z) with Q' = \J l Q\J and 
XJ l T\J C T. Latter condition is sufficient to guarantee equality U t TU = T. 
If T is given by a set of generating quadratic forms or inequalities, we can 
easily check computationally if this condition is satisfied. The same is true 
for the computation of T -automorphisms of Q, which are given by all U € 
GL d (Z) with Q = U l QU and U l TU C T. 

In contrast to the classical theory, finiteness of T-perfect forms up to T- 
equivalence may be lost (cf. [JS94J). However, although possibly not finish- 
ing in finitely many steps, we can generalize Voronoi's algorithm to a graph 
traversal search of T-equivalent T-perfect forms. Here two T-perfect forms 
are called T-contiguous if they are connected by an edge of the Ryshkov 
polyhedron 1ZP\T. 

In case of T = Tq, there exists only finitely many G-perfect forms up 
to scaling and G-equivalence due to a theorem of Jaquet-Chiffelle |Jaq95| . 
So in this case we obtain a Voronoi algorithm and have the possibility to 
enumerate (in principle) all G-perfect forms up to G-equivalence. 

In general, we can apply the procedure described in Algorithm [3] with 
respect to some given linear subspace T. If the computation finishes, we 
have a proof that there exist only finitely many T-inequivalent T-perfect 
forms. 



There are a few differences to Voronoi's Algorithm [TJ One phenomenon 
that does not occur in the classical theory is the possible existence of dead 
ends. These occur at T-perfect forms Q, whenever one of the extreme rays 
R of Vt(Q) (as in (I25h ) is positive semidefinite. In this case there is no 
T-contiguous T-perfect form on the ray {Q + aR : a > 0}. In fact, the ray 
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Input: Dimension d and a linear subspace T of S d . 

Output: A complete list of T-inequivalent T-perfect forms in S d n T. 

Start with a T-perfect form Q. 

1 . Compute Min Q and describing inequalities of polyhedral cone 
(25) Vt(Q) ={Q' £T : Q'[x] > for all x G MinQ} 

2. Enumerate extreme rays R\, . . . , R^ of the cone VxiQ) 

3. For indefinite R4, i = 1, . . . ,k, determine T-contiguous 
T-perfect forms Qi = Q + aRi 

4. Test if Qi is T-equivalent to a known form 

5. Repeat steps 1.-4. for new T-perfect forms 

Algorithm 3. Voronoi's algorithm with respect to a linear subspace T. 



is in this case contained in an unbounded face of the Ryshkov polyhedron 

n. 

Another difference to the classical algorithm is that usually we do not 
know a starting T-perfect form a priori. We can however find such a form 
starting from an initial PQF Q$ in T by applying an adapted version of 
Algorithm [5J We first compute a maximal linear subspace Lq in Vt(Qq) ( as 
in (|25p ). If it is trivial, Qo is perfect. Otherwise we choose a form R in 
Lq which is not positive semidefinite. We then can apply Algorithm [2] to 
Q = Qo and R and obtain a p > such that Qi = Qo + pR satisfies X(Qi) = 
\(Qo) and MinQo C MinQi MinQo- The maximal linear subspace L\ 
in Vt{Qi) is strictly contained in Lq. By applying this procedure at most 
dimT times, we obtain a T-perfect form Q. 

Note that our viewpoint on T-perfect forms in this article differs from 
the usual one: T-perfect and G-perfect forms are usually defined via normal 
cones o f faces of K n T in T (cf. [Mar03j . |BMS92j . |Jaq95| , |Opg95| and 
OpgOlQ. A face F of 1Z is uniquely characterized by the set 

MinT = {x G Z d : Q[x] = 1 for all Q G F}. 

The normal cone of F is the Voronoi domain cone{ xx l : x G Min F} and the 
normal cone of the face F n T in T is obtained by an orthogonal projection 
of this Voronoi domain onto T. If different inner products are used, the 
resulting cones may differ, as seen in the cases of |Jaq9 5[ and |Opg95 . 



8. Eisenstein, Gaussian and Hurwitz quaternionic perfect 

forms 

As examples for the G-theory described in the previous section, we con- 
sider three cases that have been studied intensively before. 

Eisenstein forms. If d is even, then a Q G S d is said to be an Eisenstein 
form if it is invariant with respect to a group G C GLrf(Z) of order 3 acting 



16 



ACHILL SCHURMANN 



d 


2 


4 


6 


8 


10 


dimT^ 


1 


4 


9 


16 


25 


# £-perf. forms 


1 


1 


2 


5 


1628 


maximum 5 


0.9069 . . . 


0.6168... 


0.3729... 


0.2536... 


0.0360... 



Table 3. Number and maximum densities of ^-perfect forms. 




Figure 1. Voronoi graphs for ^-perfect forms for d = 4,6, 8. 
fixed-point-free on Z rf \ {0} by z i— ► Uz. For example 

G=(id d/2 ®ft -l))' 

where (g) denotes the Kronecker product. The terminology comes from the 
fact that a corresponding lattice L C M. d can be viewed as a complex lattice 
of dimension d/2 over the Eisenstein integers 

£ = / a + fte*rt/3 : a, 6 G Z | , 

that is, L = B£ d / 2 C C d / 2 with a suitable S G GL d/2 (C). On the other hand, 
it can be seen that each complex lattice of this form yields an Eisenstein 
form. 

It turns out that the space of G-invariant forms Tq has dimension (d/2) 2 . 
In particular for d = 2 we find only one Eisenstein form up to scaling, 
associated to the hexagonal lattice A 2 . It is trivially ^-perfect (Eisenstein 
perfect). From dimension 4 on the situation is already more interesting. In 
Table [3] we list number of classes and maximum sphere packing densities of 
£ -perfect forms up to dimension 10. Figure [U shows the found contiguities 
up to dimension 8 

For d = 4, the Ryshkov polyhedron is 4-dimensional in S d (which itself has 
dimension 10). Up to ^-equivalence (by mappings Q t— > U l QU preserving 
Tq), there is only one £ -perfect form, namely the one associated to the lattice 
D4. Consequently the Voronoi graph (up to ^-equivalence) is just a single 
vertex with a loop. In dimension 6, we find already two ^-inequivalence 
^-perfect forms, associated to the lattices Eg and its dual Eg. 

The classification of Eisenstein forms in dimension 8 was almost finished 
by Sigrist in j ~Sig04| . He found all five classes of ^-perfect forms and their 
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D 4 D 6 E 8 D 8 

Figure 2. Voronoi graphs for (/-perfect forms for d = 4,6,8. 

neighboring relations. However, he could not rule out the existence of other 
£ -contiguous neighbors of the forms associated to Eg. Recently we finished 
the classification using a C++-implementation of the algorithms described 
in Sections S] and [3 The forms labeled P13, P43 and Pq2 in Figure Q] are also 
perfect forms in the classical sense. The index of the labels corresponds to 
the number of the class given in the complete list of 8-dimensional perfect 
forms that can be obtained from our webpageQ The lattice associated to Pq2 
is also known as Barnes lattice Lg (see |Mar03j Section 8.4]). The "Sigrist 
form" labeled Si is an example of an ^-perfect form which is not perfect in 
the classical sense (as already observed in |Sig04| ). 

Using our implementation we were also able to enumerate all 10-dimensional 
^-perfect forms, showing that their total number "explodes" to 1628. The 
data of our classification can be obtained from our webpageH The files 
contain a complete description of the Voronoi graph. 

Note that the largest known lattice sphere packing density 5 is attained 
among ^-perfect forms up to dimension 10. A noteworthy phenomenon that 
occurs among these forms in dimension 10 is the existence of £-inequivalent 
^-perfect forms, which are nevertheless arithmetically equivalent. This hap- 
pens for two arithmetically equivalent forms associated to the lattice K[ 
(see |Mar03| Section 8.5]). 

Gaussian forms. For even d, a Gaussian form Q E 5> is defined as a 
form containing a group G C GL^(Z) of order 4 in their automorphism group 
acting fixed-point-free on 7L d \ {0}. For example 

G =( id ^(! "o 1 ))- 

A corresponding lattice L C M d can be viewed as a complex lattice of di- 
mension d/2 over the Gaussian integers 

g = {a + bi : a,b eZ} . 

Vice versa, every such lattice yields a Gaussian form. 

As in the case of Eisenstein forms, it turns out that Tq has dimension 
(d/2) 2 . For d = 2 we find only one (/-perfect (Gaussian perfect) form up 
to scaling, namely Z 2 . As shown in Figure [21 the only (/-perfect forms in 
dimension 6 and 8 are associated to the lattices Dg, Dg and Eg. As shown in 

1 see http: //fma2 .math. uni-magdeburg . de/~achill/perfect- forms- dim8 . txt 

2 see |http: // fma2 .math .uni-magdeburg. de/~achill/E-perf ec t-f orms-dim?? .txt | 
where ?? should be replaced by 4, 6, 8 or 10. 
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n 


2 


4 


6 


8 


10 


dimTg 


1 


4 


9 


16 


25 


# (/-perf. forms 


1 


1 


1 


2 


> 17757 


maximum 5 


0.7853... 


0.6168... 


0.3229. . . 


0.2536... 





Table 4. Number and maximum densities of (/-perfect forms. 



d 


4 


8 


12 


16 


dim T-ft 


1 


6 


15 


28 


# W-perf. forms 


1 


1 


8 


7 


maximum 5 


0.6168... 


0.2536... 


0.03125... 


0.01471 . . . 



Table 5. Number and maximum densities of W-perfect forms. 

Table H] the number of equivalence classes (/-perfect forms in dimension 10 
grows even beyond the corresponding number for ^-perfect forms. So far we 
were not able to finish the classification, but we think it is computationally 
within reach on a suitable computer. 

As in the case of ^-perfect forms, the enumeration in dimension 8 was 
started by Sigrist |Sig04|. However, he did not finish the classification of Q- 
contiguous (/-perfect neighbors of Eg. Nevertheless, our computations show 
that his list was nevertheless complete. The data of our classification can 
be obtained from our webpageH Note that in dimensions not divisible by 4, 
the forms giving the densest known lattice sphere packing are not Gaussian. 

Hurwitz quaternionic forms. For d divisible by 4, a form Q E 5> is 
called Hurwitz quaternionic if it is invariant with respect to a group G C 
GL^(Z) isomorphic to 2A/± and acting fixed-point-free on Z rf \ {0}. Here, 
denotes the alternating group of degree 4. There is a correspondence 
between Hurwitz quaternionic forms and lattices in W 1 which can be viewed 
as Hurwitz quaternionic lattices over the Hurwitz quaternionic integers 

H = {a + bi + cj + dk : a,b,c,d E Z or a,b,c,d E Z + 5} . 

We refer to [CS991 Section 2.6] for details. 

It turns out that Tq is of dimension { d (^)- This leaves only one Hurwitz 
quaternionic form and therefore only one TC -perfect form up to scaling for 
d = 4, which is associated to D4. As shown in Figure there is also only 
one equivalence class of W-perfect forms in dimension 8, corresponding to 
Eg. 

The situation becomes more interesting in dimension 12 (cf. Table [SJ. 
By our computations, there are precisely eight classes of W-perfect forms, 
as previously observed by Jaquet-Chiffelle and Sigrist (cf. |Sig08|). Figure [3] 
uses their labeling. The data of our computations can be obtained from 

^ see |http : //fma2 .math.uni-m agdeburg. de/~achill/G-per f ect-f orms-dim?? . txt | 
where ?? should be replaced by 4, 6 or 8. 
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Qmax 




MAX 



u 

Figure 3. Voronoi graphs for "H-perfect forms for d = 8, 12. 

our webpageEl Note that all W-perfect forms are also perfect in the classical 
sense. 

A quite interesting consequence of the classification in dimension 12 is the 
possibility to derive of a sharp bound for the largest possible sphere packing 
density among Hurwitz quaternionic forms in dimension 16, as shown by 
Vance [Van09j using a Mordell type inequality. She shows that the Barnes- 
Wall lattice BWi6 has the largest density among lattices with a Hurwitz 
quaternionic structure in dimension 16. A very nice example of a human- 
computer-interacted proof! 
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GAP — Groups, Algorithms, Programming - a system for computational dis- 
crete algebra, ver. 4.4., http://www.gap-system.org/ 
Irs by D. Avis, ver. 4.2., http : //cgm. cs .mcgill . ca/~avis/C/lrs .html 
MAGMA — high performance software for Algebra, Number Theory, and Geom- 
etry, ver. 2.13., http://magma.maths.usyd.edu.au/ 

shvec by F. Vallentin, ver . 1.0., see http: //www.math.uni-magdeburg . de/lattice_geometry/ 
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